Local-channel-induced rise of quantum correlations in continuous-variable systems 
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It was recently discovered that the quantum correlations of a pair of disentangled qubits, as measured by the 
quantum discord, can increase solely because of their interaction with a local dissipative bath. Here, we show 
that a similar phenomenon can occur in continuous-variable bipartite systems. To this aim, we consider a class 
of two-mode squeezed thermal states and study the behavior of Gaussian quantum discord under various local 
Markovian non-unitary channels. While these in general cause a monotonic drop of quantum correlations, an 
initial rise can take place with a thermal-noise channel. 
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I. INTRODUCTION 

Within the remit of quantum information processing (QIP) 
theory fUl and beyond, a paramount topic is the study of 
quantum correlations (QCs). Until recently, the emergence of 
QCs has been regularly highlighted in connection with non- 
separable states of multipartite quantum systems j^. The vi- 
olation of the celebrated Bell inequalities is a typical signa- 
ture of the extra amount of correlations that quantum systems 
can possess besides those of a purely classical nature 01 [H]. 
Non-separability, i.e. entanglement, which relies on the su- 
perposition principle, has been in fact regarded as a neces- 
sary prerequisite in order for QCs to occur. In 2001, however, 
it was discovered |3j] a new way in which the superposition 
principle can entail an alternative type of non-classical corre- 
lations even in the absence of any entanglement. Typical in- 
stances are mixtures of pure separable states which are locally 
non-orthogonal, i.e. indistinguishable U More rigorously, 
a state has non-classical correlations whenever its associated 
density operator cannot be diagonalized in a basis which is 
the tensor product of local orthonormal bases J5|. When this 
occurs the entire correlations' content cannot be retrieved by 
any local measurement at variance with the classical frame- 
work where this is always achievable. Following such break- 
through, a growing interest has developed especially after the 
realization that this new QCs' paradigm could be key to some 
known entanglement-free QIP schemes ^ . 

Various measures have been proposed in the literature to 
detect such kind of QCs. A prominent one is quantum discord 
whose definition simply arises from the quantum general- 
ization of two classically-equivalent versions of the mutual in- 
formation (one being based on the conditional entropy). These 
are found to differ in the quantum framework, the correspond- 
ing discrepancy just being the quantum discord. This encom- 
passes both QCs associated with entanglement and those that 
can be exhibited by separable states. The non-equivalence 
between such two forms of non-classical correlations arises 
when mixed states are addressed since as long as the state is 
pure the discord coincides with the entropy of entanglement. 

Significant motivations to pursue a deeper understanding 
of these problems are currently being provided by various in- 
vestigations targeting non-unitary dynamics. While entangle- 
ment is known to be extremely fragile to environmental inter- 
actions, QCs - as given by such extended notion - are gen- 
erally very robust and in some cases even fully insensitive to 



phase noise 01 • This is related to the fact that zero-discord 
states are a set of negligible measure within the entire Hilbert 
space [8]. Still in the framework of non-unitary dynamics, 
very recently it was found that the interaction with a local and 
memoryless bath can even create QCs initially fully absent 
!9- [l2ll . Such effect, which is unattainable with entanglement 
1 2[], was shown to take place in particular for qubits prepared 
in a fully classical state and undergoing a local amplitude- 
damping channel. This describes the dissipative interaction 
with a local bath (such as the spontaneous emission of a two- 
level atom) [lj. Specifically, the behavior consists of an ini- 
tial rise of QCs until a maximum is reached followed by a 
slow decay, the entanglement being zero throughout. The es- 
sential underlying mechanism is that such non-unitary process 
can map orthogonal states of a subsystem onto non-orthogonal 
ones. The state thus can no more be diagonalized in a tensor 
product of orthonormal bases, i.e. it acquires QCs l9l [TTll . 

Most of the work carried out so far along these lines, 
though, tackled qubits. It is however natural to wonder how 
the above findings are generalized for quantum continuous- 
variable (CV) systems, which routinely occur in quantum op- 
tics. Such investigations are still in their infancy. Indeed, the 
first measure of generalized QCs for Gaussian states of bipar- 
tite CV systems has been worked out only last year lfl3l [Till . 
Such quantity, called Gaussian discord, is basically defined 
in the same spirit of standard discord [ 3fl and likewise can 
be non-zero for separable states lfl3l [l4 l. Recently, another 
measure has been proposed lfl5ll . To date, only few studies 
lfl3i [l6tl have targeted the non-unitary dynamics of Gaussian 
discord in the presence of environmental interactions includ- 
ing one model of non-Markovian reservoir lfl8ll . It was found 
that apart from the case of a common reservoir the non-unitary 
dynamics is detrimental to QCs II13[ Il6l Il7ll . 

Our goal in this work is to analyze the behavior of Gaus- 
sian discord in the case of two CV modes under the most rel- 
evant known instances of local - i.e. single-mode - Gaussian 
memoryless channels lTl9l - l2lll each described by an associated 
completely positive quantum map. Our major motivation is to 
assess whether growth of QCs can take place for some local 
non-unitary channels as in the case of qubits lr9l-HTll. We will 
indeed find that in significant analogy with the qubits' frame- 
work this can occur with the lossy channel (more in general 
with low-temperature thermal-noise channels). 

This paper is organized as follows. In Section[n] we briefly 
review the definition of Gaussian discord. In Section [HI] we 
describe the class of squeezed thermal states, which is the one 
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which will be our focus in this work. In Section|IV] we review 
each of the aforementioned Gaussian channels and the related 
main features that we will refer to. In Section [V] we show 
the behavior of QCs under each of the considered channel for 
some paradigmatic initial states. In Section [VlJ we provide 
a simple picture that allows to obtain insight into the effects 
presented in the previous section. Finally, in Section [VIII we 
draw our conclusions. 



II. QUANTUM DISCORD OF GAUSSIAN STATES: 
REVIEW 

Given two systems 1 and 2 in a state p, the quantum mutual 
information is defined as 



I(p)=S(p 1 )+S(p 2 )-S(p) 



(1) 



where pip) = Tr2(ijp is the reduced density operator describ- 
ing the state of 1 (2) and S(cr) = -Tr(crlogcr) is the Von 
Neumann entropy of an arbitrary state <x. A local general- 
ized measurement on 2 can be specified by a complete set 
of positive-operator-valued projectors (POVM) {IT}, where 
k indexes a possible outcome. If k is recorded with prob- 
ability pi< = Tr[pl®n<.] the overall system collapses onto 
the (normalized) state pk = (pi ®llk)lpk- The maximum of 
S(Pi) — Yik PkS(Pk), i- e - the mismatch between the entropy of 
1 and the average conditional entropy, reads 



J"(p) = max 

{n 4 | 



S(pi)-^PkS(pk) 



(2) 



and is taken over all possible POVM measurements each de- 
scribed by the operator 11^. If 1 and 2 are CV modes and 
one restricts to generalized Gaussian measurements ll22Tl the 
Gaussian discord D*~ is defined as the discrepancy between 

©and© mm 



IT 



I- J 



(3) 



The arrow reminds that measurements over mode 2 have been 
considered. 

As is well-known, any two-mode zero-mean Gaussian 
state is fully specified by its covariance matrix cry = 

Tr[pi2 (RiRj+RjR;)], where R = {x\,p\,x 2 ,p 2 \ is the set of 
operators corresponding to the phase-space coordinates. Up 
to local symplectic, i.e. unitary, operations the covariance ma- 
trix can be arranged in the form 



A C 
C B 



(4) 



where A = diag(a, a), B = diag(£>, b) and C = diag(ci , C2). The 
determinants of such diagonal matrices I\ = det A, I 2 = det B, 
I3 = detC along with I4 = deter are symplectic invariants, 
i.e. they are invariant under local symplectic transformations. 
The quantum discord is a function of these symplectic 
invariants according to lfl3l [Till 



D^=h(ylh)- ll(d-) -h(d + ) + hi^ 



a/77 + 2 ^Ihh + 2/3 

1 + 2 vt; 



(5) 



where h(x) = (x+ 1 /2) ln(x+ 1 /2) - (x- 1 /2) ln(*- 1 /2) and 



d 2 _ A ± VA 2 - 4/4 



(6) 



with A = 1 1 + h + 21 j. The discord in terms of measurements 
on mode 1 D"* is obtained from © upon the replacement 
I\ <-> h- In the remainder of this work, we will be solely 
interested in the discord ©, hence we will drop the subscript 
henceforth and set D = . 



m. TWO-MODE SQUEEZED-THERMAL STATES 

In this work, we shall focus on the class of two-mode 
squeezed-thermal states (STSs) whose generic element reads 



r(a ]a*-a 1 a 2 ) 



(Pi ®P2) 



(7) 



p = e 

where each single-mode thermal state p,- (i = 1, 2) is given by 

N? 

rT 

n=0 



Pi= / r\n)i(n\ 

H ^ (1 +Nd" +1 



(8) 



with Nj the corresponding average number of photons. 

The characteristic function x(M, ^2) lf23ll is obtained from 
© as 



X(A u A 2 )=TT\pD(A 1 )D(A 2 )] , 



(9) 



where D{Aj) = exp(/l,fl / ! - /Pa,) is the displacement oper- 
ator of the /th mode (; = 1,2 with a, and a] the usual 
ladder operators of the /th mode). The density opera- 
tor is retrieved from the characteristic function as p = 
/ d 1 Aid 2 A 2 x{AuA 1 )D(-Ai)D( r A 2 )ln 1 . It is then straightfor- 
wardly checked that 

^(/ll Ao)~ cosh/-/li-sinh/-/l,| 2 ^-(A^2+l/2)| cosh rArsinh rA\ | 2 

(10) 

As states © are Gaussian, the covariance matrix is related to 
the characteristic function according to l24ll 



X(A) = exp 



A T crA 



(11) 



where A=(ori,/?i, a2,/3z) an d we have decomposed each 
complex variable Aj (j = 1,2) as Aj = (aj + iBj)/ y2. Upon 
comparison between Eqs. ( [Tol l and (fTTT > and in the light of ©, 
the diagonal matrix elements defining A, B and C [cf. Eq. ©] 
are obtained as 

a = (1 +N r )Ni +N r N 2 + N r + 1 /2 , (12) 
b = N r N\ + ( 1 +N r )N 2 +N r +l/2 , (13) 
c l = -c 2 = -(l+Nl+N 2 )^fNAT+Nr) , (14) 

where we have set N r = (sinh r) 2 . 
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IV. LOCAL GAUSSIAN CHANNELS 

In this Section, we review the salient features of the lo- 
cal one-mode Gaussian channels whose effect on the quan- 
tum correlations of two-mode states we aim to scrutinize. 
Throughout, we will assume that each of such channels acts 
on mode 2 only. The channels are Gaussian in that each of 
these maps Gaussian states into Gaussian states, hence allow- 
ing for use of Gaussian discord lfl3l [Till to shed light onto 
the QCs' behavior. Specifically, we address the thermal-noise 
channel, which reduces to the lossy channel in the limit of 
vanishing reservoir temperature, the amplifier channel and the 
classical-noise channel rf20ll . 



where k > 1 is a gain parameter. In the limit k — 1 the channel 
reduces to the identity operator. By proceeding analogously 
to the case of the thermal-noise channel we find that the co- 
variance matrix has again the same form as in Eq. <j4j with the 
diagonal matrix elements now given by 

a! = a , (20) 
b' = kb + (k-l)(N+^\ , (21) 

c', = -c' 2 = Vici . (22) 
C. Classical-noise channel 



A. Thermal-noise channel 

This channel describes the dissipative interaction of a 
single-mode CV system with an environment at thermal equi- 
librium. Indeed, its associated map is routinely worked out by 
assuming that a single-mode environment in a thermal state 
specified by the average photon number N is mixed by a beam 
splitter with the single-mode system. The state of the lat- 
ter, and hence the map describing the channel, can then be 
obtained by simply tracing out the environmental degree of 
freedom. The channel quantum efficiency is measured by the 
parameter 77 such that < 77 < 1 (in the above model 77 is the 
transmissivity of the beam splitter). 

Under the thermal-noise channel of efficiency 77 and envi- 
ronment's average photon number N, the characteristic func- 
tion x(A\, At) transforms into^'(/li, A2) according to i25ll 

xxXAuA 2 )= X {M,^A 2 )e- ( - l -^ N+ ™? (0<i;<1). (15) 

For N-Q (zero-temperature environment) Eq. ( fTBI l reduces to 
the case of a lossy channel I20I1 acting on mode 2. 

By substituting ( TTOb in Eq. ( fTTT i and comparing this with 
Eq. (fT5T l it is straightforwardly found that the covariance ma- 
trix keeps the same structure as in Eq. (0]i with A = diag(a', a'), 
B = diag(/b', b') and C = diag(c'j , c' 2 ). The new diagonal matrix 
elements {a', b', c'} are related to the input ones [cf. Eqs. (fT2l - 
(fT~4T >1 according to 

a! = a , (16) 

b' =nb + (l-r])^N+^\ , (17) 

c\ = -c' 2 = V^C! . (18) 

B. Amplifier channel 

This channel shares features similar to the thermal-noise 
channel but with the essential difference that it brings about an 
intensity amplification instead of a damping. Its correspond- 
ing map changes the characteristic function according to l26ll 

X , {AuA 2 )=x{Au^kA 2 )e< k -^™^ (k>l) . (19) 



This channel arises when classical Gaussian noise is super- 
imposed on the single-mode system The parameter on 
which it depends is the number of injected noise photons n > 
in a way that the corresponding map becomes the identity in 
the limit « = 0. The characteristic function is transformed un- 
der this channel according to 

x'(M,A 2 )= X (A u A 2 )e-"^ 2 (*>0). (23) 

By comparing this with Eqs. (TTOb and ( fTTT i it turns out that 
the covariance matrix has the form (01 with diagonal matrix 
elements 

a' = a , (24) 
b' = b + n , (25) 
c\ = -c' 2 = ci . (26) 

V. BEHAVIOR OF QUANTUM CORRELATIONS 

Our aim in this Section is to present some typical behav- 
iors of Gaussian discord when states of the family (0 are sub- 
ject to the local Gaussian channels introduced in the previous 
Section. Unfortunately, a general analysis is quite demanding 
owing to the complicated functional form of (0. Nonethe- 
less, as will become clear, one can obtain significant insight 
by addressing some paradigmatic instances. To begin with, 
we consider initial STSs [cf. Eqs. (01 and (©] having r — 1 
and iVi = N 2 . In Fig. 1, we plot the behavior of Gaussian dis- 
cord (|5]l for increasing values of iVi = N 2 in the presence of 
the thermal-noise and amplifier channel (both with = 0) as 
well as the classical-noise channel. We recall that each lo- 
cal non-unitary map affects only mode 2. The classical-noise 
channel [see Fig. 1(c)] has a merely detrimental effect on the 
QCs, which monotonically decay for a growing number of 
noise photons regardless of N\ . A similar effect is exhibited 
in the case of the amplifier channel when the gain parameter 
k is increased [see Fig. 1(b)]. In contrast, however, a non- 
monotonic behavior can take place with the lossy channel, as 
shown in Fig. 1(a). For low average photon numbers of each 
mode N\ = N 2 , similarly to Figs. 1(b) and (c), a monotonic 
drop of QCs occurs as the quantum efficiency rj decreases (we 
recall that the corresponding map reduces to the identity oper- 
ator for 77 — 1). As the photon number becomes larger, though, 
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FIG. 1. (Color online) Behavior of Gaussian discord D under a thermal-noise (a), an amplifier (b) and a classical-noise channel (c) with the 
former two taken for N = 0. Throughout, we have considered an initial state Q for r = 1 and N] = N 2 = 1 (solid black lines), N[ = N 2 = 10 
(red dashed) and Nj = N 2 = 100 (blue dotted). The inset in (a) shows D against ;/ for N[=N 2 = 1000 under the zero-temperature thermal-noise 
channel. 





FIG. 2. (Color online) Gaussian discord D against 77 under the 
thermal-noise channel for the average photon number of the reser- 
voir N = (lossy channel, solid black line), N = 1 (red dashed), 
N= 10 (blue dotted) and N = 50 (green dot-dashed). Throughout, the 
considered initial state is with r= 1 and Ni=N 2 = 10. 



FIG. 3. (Color online) Initial slope of the Gaussian discord p, 
i.e. dD/drj at 77 = 0, vs. N and N t = N 2 rj for r = 1 under the thermal- 
noise channel. The intersection line between the function profile (in 
blue) and the plane p = (in brown) specifies, for set N t , the maxi- 
mum value of N yielding QCs' rise. 



the discord remarkably undergoes a rise so as to reach a max- 
imum value and eventually drop to zero. Such increase can 
be quite significant. For instance, as shown in Fig. 1(a), in 
the case N\ = N2 = 10 the maximum taken by D is about 2.5 
times the initial value (for 77 = 1) and even ^10 larger when 
Ni=Nt= 50. For higher photons numbers (still keeping fixed 
the set value of r) the initial discord becomes extremely low 
but a significant rise still takes place before the asymptotic de- 
cay. We illustrate the latter feature in the inset of Fig. 1(a) for 
the paradigmatic cas e N\ =Nz = 1000, which corresponds to a 
fully separable state 12811 with an initial discord D\ n =\ — 10~ 6 . 
Yet, under the lossy channel this undergoes a rise reaching 
a maximum 0.5 before the eventual slow decay. For all 
practical purposes, one can thus regard this behavior as mere 
creation of (previously absent) quantum correlations, which is 
significantly reminiscent of an analogous effect occurring for 
qubits under local amplitude-damping channels lr9l- [Tlll . 

It is natural to wonder how the discord rise is affected in 
the presence of a reservoir having non-zero temperature. To 
investigate this, in Fig. 2 we display T> against 77 in the case 
of the thermal-noise channel for growing values of the reser- 
voir average photon number N. Clearly, the effect of temper- 



ature is to spoil the discord increase. At high enough temper- 
atures, the initial rise is no more exhibited and the behavior 
reduces to a mere monotonic decay similarly to the amplifier 
and classical-noise channels [cf. Figs. 1(b) and (c)]. Note that 
there is a threshold value for N separating the non-monotonic 
regime (featuring QCs rise) from the monotonic one, the lat- 
ter occurring above the threshold. For set r, such critical value 
linearly grows with N\ . This is shown in Fig. 3, where we plot 
the initial slope of the Gaussian discord p, i.e. p = dD/drj at 
77 = 0, vs. and N\ = A^- By recalling that 77 decreases dur- 
ing the system's evolution under the thermal-noise channel, 
p > (p < 0) means that an initial decay (rise) of QCs oc- 
curs. Hence, the intersection of p(N,Ni) with with the N~N\ 
plane provides the functional dependance of the aforemen- 
tioned threshold value on N\ . 



VI. INSIGHT INTO THE RISE OF QUANTUM 
CORRELATIONS 

Here, we give a picture that illustrates the emergence of 
some key features and behaviors presented in the previous 
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FIG. 4. (Color online) (a) Contour plot of Gaussian discord T) against b' and d . (b) Oriented trajectories in the b' — c' plane arising with the 
lossy channel (solid line), the zero-temperature amplifier channel (dashed) and the classical-noise channel (dotted). The red point corresponds 
to the initial state, (c) Zoom of the previous figure highlighting the initial rise of discord occurring only under the lossy channel, (d) Effect of 
temperature in the case of the thermal-noise channel: oriented trajectories for ^ = (lossy channel, solid line), N=l (dashed), N= 10 (dotted) 
and ^ = 50 (dot-dashed). Throughout, we have considered an initial state (|7) for r = 1 and Ni =N2 = 10. 



Section. We focus on the initial state specified by r — 1 and 
N\ -N2-IO exhibiting a rise of discord followed by a decay 
in the case of the thermal-noise channel (at low temperatures) 
and a monotonic drop with the other two channels (cf. Figs. 1 
and 2). As in each of such three cases the environment does 
not affect parameter a' [see Eqs. ( fTo*b . d20t and (I24H and the 
equality c\ - -cL always holds, in the addressed regime the 
discord is in fact a function of V andc' = |c'J — \c' 2 \ [cf. Eq. ((5])]. 
As shown in the contour plot in Fig. 4(a), in its domain of def- 
inition D decreases with b' and increases with c' . One can 
now obtain insight into the behavior of QCs by observing that 
the equations for b' and c' define a characteristic curve (as- 
sociated with the specific channel) parametrized by 776 [0, 1], 
k > 1 and n > in the case of the thermal-noise, amplifier and 
classical-noise channel, respectively. By expressing each of 
such parameters as a function of c' and replacing it into the 
equation for b' , one obtains the trajectory of each channel in 
the b'—c' plane as 



,, b-(N+l/2) , 2 



c' z + (N+l/2) (c' G [0,c]) , (27) 



, b+(N+l/2) , 2 



b' = 



c' z -(N+l/2) (c e [c,oo]) , (28) 



c' = c (b'e[b,oo]) 



(29) 



In Figs. 4(b) and (c), we report the above trajectories along 
with the same contour plot of discord as in Fig. 4(a). Evi- 
dently, due to the above discussed functional shape of T> a 
monotonic decrease of QCs necessarily takes place under the 
amplifier and classical-noise channels. It is also clear from 
the bottom left portion of Fig. 4(b) that for the lossy channel 
as 77 is decreased T> must eventually drop. However, the tra- 
jectory is such that over the first stage (r\ slightly below 1) the 
QCs grow [this is best illustrated by the zoom presented in 
Fig. 4(c)]. 

The effect of temperature in the case of the thermal-noise 
channel (cf. Fig. 2) can be understood by scrutinizing Eq. (|27| | 
and Fig. 4(d). For r\ — > c' tends to zero and thereby b' — > 
(N+l/2). Hence, the higher N (i.e. the temperature) the larger 
the asymptotic value of V . Accordingly, the concavity C - 



b — (N+l/2) /c 2 - which is positive for N = - progressively 
decreases so as to eventually become negative. As soon as it 
becomes small enough, the shape of D(b',c') [cf. Fig. 4(a)] 
prevents from any initial rise to take place, which results in a 
monotonic decay (see Fig. 2). 



VII. CONCLUSIONS 

In this paper, we have investigated the behavior of Gaussian 
discord for a two-mode squeezed thermal state subject to var- 
ious local Gaussian channels. Mainly motivated by the find- 
ing that a local amplitude-damping channel can create QCs 
in a pair of qubits, we have explored the dynamics under the 
thermal-noise, amplifier and classical-noise channel. While 
the typical behavior is a monotonic decrease of QCs, as one 
would expect, we have found that in significant analogy with 
the qubits' framework a thermal-noise channel can give rise 
to a non-monotonic behavior comprising an initial rise of dis- 
cord. For large enough photon numbers of each mode such 
that the initial QCs are in fact negligible, the above entails 
creation of previously absent QCs for all practical purposes. 
The reservoir temperature spoils this phenomenon in a way 
that when it is high enough a mere monotonic decay occurs. 

We have provided a picture that allows to shed light on var- 
ious features, in particular the reason why the discord rise is 
exhibited only under the thermal-noise channel as well as the 
detrimental effect of temperature on it. 

These findings significantly extend to the CV-variable sce- 
nario one of the most counter-intuitive effects entailed by the 
emerging extended paradigm of QCs: quantum correlations 
can be established in a composite system through the interac- 
tion with a local and memoryless bath. 
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